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The radioactive radium-225 (^^®Ra) atom is a favorable case to search for a permanent electric 
dipole moment (EDM). Due to its strong nuclear octupole deformation and large atomic mass, 
^^®Ra is particularly sensitive to interactions in the nuclear medium that violate both time-reversal 
symmetry and parity. We have developed a cold-atom technique to study the spin precession of ^^®Ra 
atoms held in an optical dipole trap, and demonstrated the principle of this method by completing 
the first measurement of its atomic EDM, reaching an upper limit of |d(^^®Ra)| < 5.0 x 10“^^ e-cm 
(95% confidence). 

PACS numbers: ll.30.Er, 24.80.-by, 32.10.Dk, 37.10.Gh 


The existence of a permanent electric dipole moment 
(EDM) would violate both time-reversal (T) and parity 
(P) symmetries [T], and, by the CPT theorem, charge- 
parity (CP) symmetry [2]. Although the Standard Model 
(SM) exhibits CP violation [3H5], its effects contribute to 
EDMs only at higher orders [5] , resulting in EDM values 
too small to be detected in the foreseeable future. On the 
other hand, extensions of the SM, such as supersymme¬ 
try, naturally provide additional sources of CP violation 
that explain the matter-antimatter asymmetry of the uni¬ 
verse and result in observable EDMs [S]. Limits derived 
from recent EDM searches in the neutron the i^^Hg 
atom [5], and the electron (ThO) [5] provide some of the 
most stringent limits for these beyond-SM CP-violating 
interactions Ha El- 

Searches in diamagnetic atoms (^®®Hg [8], ^^®Xe m) 
and molecules (TIE [T^) are primarily sensitive to CP- 
violating interactions in the nucleus M- Although the 
nuclear EDM is shielded by the atomic electrons, this 
shielding is imperfect due to the finite size of the nu¬ 
cleus |15j . The remaining measurable quantity, charac¬ 
terized by a radially-weighted EDM of the nucleus, is 
called the nuclear Schiff moment, which induces a pro¬ 
portional and measurable atomic EDM. The best limits 
on CP-violating interactions originating within the nu¬ 
cleus are derived from the limit on the atomic EDM of 
i99Hg: |d(i99Hg)|< 3.1 X 10“^® e-cm [16]. Here, we re¬ 
port the first measurement of the EDM of the diamag¬ 
netic ®^®Ra atom. P-odd and T-odd effects in ^®®Ra are 
greatly enhanced by collective effects and a closely spaced 
parity doublet structure resulting from the octupole de¬ 
formation of its nucleus m- Indeed, nuclear structure 
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calculations [Ill[l9], combined with atomic structure cal¬ 
culations HO], predict the EDM of ^®^Ra to be 2-3 orders 
of magnitude larger than that of i®®Hg E]. This en¬ 
hancement factor may be even larger given the sizeable 
variations in the calculations of ^®®Hg m- Meanwhile, 
these special nuclear structure properties also allow the 
Schiff moment of ^^®Ra to be more reliably calculated 
than that of i®®Hg [TO] . 

®^®Ra is also attractive from an experimental perspec¬ 
tive. The 14.9 day half-life allows this isotope to be 
obtained as a radioactive source in sufficient quantities 
for experiments to run off-line, away from an accelera¬ 
tor. Its nuclear ground state has spin I = 1/2, which 
eliminates some systematic effects and decoherence aris¬ 
ing from electric field gradients. Finally, the atomic 
structure of radium enables laser cooling and trapping 
[5^ [53]. To search for an EDM, spin precession frequen¬ 
cies (w±) under the influence of both a uniform electric 
field (E) and a uniform magnetic field (B) are measured: 
huj± = 2/rP ± 2dE. Here d is the EDM; /r is the magnetic 
dipole moment (— 0.7338(15)/rAr |21|); are precession 
frequencies corresponding to the E-field parallel or anti¬ 
parallel to the B-field. The EDM measurements were 
performed with ^®®Ra atoms in a standing-wave optical 
dipole trap (ODT) (FigQ. This approach was theoreti¬ 
cally examined in detail for ^®®Hg [25] . 

®^®Ra is a decay product of the long-lived ^^®Th iso¬ 
tope {ti /2 = 7300 yr), and can be chemically separated 
from a ^^®Th stock sample. The U.S. National Isotope 
Development Center can provide 10-12 mCi of ^®^Ra ev¬ 
ery two months |26| . For the two experimental runs pre¬ 
sented in this paper, ^®^Ra samples of 3 mCi and 6 mCi 
(10^“^ atoms), respectively, were loaded into an oven, and 
used gradually over the following two weeks. In addition, 
the long-lived, naturally-occuring ^®®Ra (ti /2 = 1600 yr) 
was a useful proxy during the development phase of the 
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experiment, and for tuning the apparatus prior to the 
EDM measurements. A typical load of ^^®Ra is 2-5 /rCi, 
amounting to 10^® atoms, or 100 times more atoms than 
the typical ^^®Ra sample. However, with 1 = 0, ^^®Ra 
cannot be used in a nuclear spin precession measurement. 

Laser manipulation of radium atoms has previously 
been described in detail [53]. Neutral radium atoms, 
both ^^®Ra and ^^®Ra, pass through a transverse cool¬ 
ing region and a Zeeman slower, before being captured 
in a three-dimensional (3D) magneto-optical trap (MOT) 
based on the intercombination transition 7s^ ^Sq F = 1/2 
—>■ 7s7p ^Pi F = 3/2 at 714 nm. With an associated life¬ 
time of 422 ± 20 ns this transition is weak compared 
to those typically used for primary cooling and trapping 
of atoms. However, it is the transition in radium that 
is the nearest to being closed, with a branching ratio of 
only 4x10“^ leaking to 7s6d ^Di [28l|29]. Atoms in 7s6d 
^Di are re-pumped to the ground state with a 1429 nm 
laser. Largely due to the weakness of the laser force, the 
trapping efficiency, from the sample in the oven to the 
MOT, is only 1x10“®. Typically, 10® ^^®Ra atoms or 
10® ^^®Ra atoms are accumulated over the MOT lifetime 
of 40 seconds and cooled to 40 ± 15 ^K. 

Next, the atoms are transferred to an optical dipole 
trap (ODT) formed by focusing a 40 W, 1550 nm, hori¬ 
zontal laser beam with a 2 m focal length lens. The re¬ 
sulting trap is elongated: 1 cm long and 100 ^m across, 
with a depth of 500 /iK. The transfer efficiency from the 
MOT to the ODT is typically 80%, helped by the fact 
that 1550 nm is nearly magic to the cooling transition, 
meaning the differential light shift of the ground and ex¬ 
cited states is comparable to the natural linewidth [55] . 
The ODT is then translated axially over a distance of 1 m 
into a separate chamber. Here it overlaps a perpendicu¬ 
lar standing-wave ODT formed by a single-mode, single¬ 
frequency, 10 W, 1550 nm retroreflected laser beam fo¬ 
cused down to 100 /im in diameter, and linearly polar¬ 
ized in the horizontal direction. The transfer of atoms 
from the first ODT to the second is assisted by briefly 
turning on a one-dimensional (ID) MOT to compress 
the aspect ratio of the atom cloud from 100:1 to 1:1. 
This 60 fim cloud retains its shape in the standing-wave 
ODT after the first ODT is turned off. The overall ef¬ 
ficiency of transferring atoms from the 3D MOT to the 
standing-wave ODT is 5%. No residual magnetization of 
the shields (see below) was detected after the ID MOT 
was switched off; any residual magnetic field was found 
to be less than the measurement sensitivity of 300 nG. 

Figure [2 provides the layout surrounding the atoms in 
the standing-wave ODT. The trap is placed at the center 
of a pair of parallel copper electrodes, cylindrical with 
a vertical axis, whose end faces are 1.6 cm in diameter 
and 2.3(1) mm apart. The upper electrode is grounded; 
the lower one can be ramped to voltages between -1-15.5 
kV and —15.5 kV, generating a uniform E-field of 67 
kV/cm in either the up or down direction (parallel or 
anti-parallel to the B-field, respectively). The leakage 
current measured on the grounded side is typically < 
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FIG. 1. (Color online) Diagram of the center of the science 
chamber. A standing-wave ODT is aligned between a pair 
of copper electrodes. A collimated beam at 483 nm both 
optically pumps and probes the atoms, providing absorption 
images of the atom cloud. A traveling-wave ODT propagating 
in the direction normal to the page delivers atoms from the 
3D MOT to the standing-wave ODT. Inset: CCD image of 
atoms loaded in the standing-wave ODT; the image is 450 pm 
by 450 pm. 


80 pA. The variation of the E-field near the center is 
<l%/mm. The electrode assembly is inside a glass, non¬ 
magnetic vacuum enclosure, which in turn is surrounded 
by a cosine-theta coil wound on an aluminum cylinder 
of 0.32 m diameter and 0.65 m length. Three layers of 
mu-metal shields surround this assembly and reduce the 
influence of external B-fields by a factor of 2 x 10'^ when 
measured in the vertical direction at the center. The coil 
inside the shields generates a stable and uniform B-field 
of 15-30 mG in the vertical direction (intentionally var¬ 
ied between experimental runs). Its spatial variation is 
<l%/cm; its instability is <0.01% when averaged over a 
load-measurement cycle of 50 s. 

The number of atoms in the standing-wave ODT is 
probed by a blue laser beam tuned to the resonance of 
7 s 2 iSo F = 1/2 ^ 7s7p ®Pi F = 3/2 at 483 nm, co¬ 
propagating along the standing-wave ODT laser beams. 
This transition can cycle for an average of about 1000 
times before leaking to lower D states. The shadow im¬ 
age produced by resonant absorption is cast onto a GGD 
camera (Fig. [H. To produce an absorption image (given 
about 200 ^^ma atoms in the trap), the blue beam is 
pulsed for 1.45 ms, during which each atom absorbs on 
average 100 photons. A linear combination of images 
without atoms is used for background subtraction; i.e., 
to suppress distortions arising from interference effects 
[30] . The noise of this detection scheme is approximately 
1.2 times the photon number shot noise. This detection 
method is destructive in that it heats the atoms out of 
the trap and pumps them to metastable levels. By mea¬ 
suring the number of atoms at various delay times, the 
lifetime of the atoms in the trap is determined to be 3-5 
s, about 10 times shorter than that in the 3D MOT. This 
is consistent with the higher vacuum pressure observed 
in the glass chamber. 
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FIG. 2. (Color online) Pulse sequence used in the EDM mea¬ 
surement. Half-period denotes half of the spin precession pe¬ 
riod. 


The spin polarization is both produced and detected 
via optical pumping by a circularly polarized blue beam 
tuned to the resonance of the non-cycling transition 7s^ 
^So F = 1/2 —>■ 7s7p ^Pi F = 1/2. An atom with a spin in 
the fully polarized state does not absorb photons, while 
one in the opposite-spin state absorbs on average three 
photons before becoming fully polarized. Compared to 
the atom-number measurement on a cycling transition, 
the spin-sensitive detection scatters 30 times fewer pho¬ 
tons, resulting in a reduction of image contrast. 

The pulse sequence used for the EDM measurement 
is shown in Fig. Two kinds of pulses, generated by 
acousto-optical modulator (AOM) switches, are used: 1.5 
ms pulses to polarize the atom cloud, and 60 /xs pulses 
to measure the number of atoms in the opposite-spin 
state (short pulses optimize the image signal-to-noise ra¬ 
tio, whereas longer ones for optical pumping maximize 
atom polarization). After each pulse, the 483 nm laser is 
blocked within 1 ms by an additional mechanical shutter 
to prevent decoherence induced by light leaking through 
the AOM while the atoms process. Four images of atoms 
are recorded in each measurement cycle. The first oc¬ 
curs prior to any polarization pulse and, thus, produces 
a signal at half of the maximum contrast. The second 
occurs after half of the spin precession period, following 
a polarization pulse, and has a signal at the maximum 
contrast. The third is taken about 2 s after polarization; 
it is during this time that the E-field is applied. Then, 
the atoms are repolarized, followed by the fourth image 
also occuring about 2 s after polarization, but this time 
with no applied E-field. The third and fourth images are 
normalized to the second one in order to reduce sensitiv¬ 
ity to atom number fluctuations. The third image is used 
to build the “E-held on” spin precession curve, and the 
fourth builds the “E-field off’ spin precession curve. The 
data indicate that the decreasing contrast is consistent 
with the lifetime of the atoms in the ODT. 

Figure [^presents the spin detection data with the E- 
field on, pointing parallel to the B-field, anti-parallel, or 
with the E-field off as functions of the free precession 
time. The E-field is at its maximum value of 67 kV/cm 
for 1.2 s during the 2 s of precession. Each data point was 
corrected for trap losses using the measured trap lifetime. 
Data taken under the three E-field configurations were 


simultaneously ht to a combined set of equations: 

2/E-fieldOff = ~ PcOs{u!t)] 

2/Parallel, Anti—Parallel — ^ p p PcOs(j^t -j- d it A0/2)] 

Five parameters A (normalization), P (atom polariza¬ 
tion), oj (precession frequency with E-field off), 9 and 
(defined below) were fit without constraints. ?/E-fieidOff 
is the integrated signal in Image ^4 after normaliza¬ 
tion to Image ^2 and background subtraction. Similarly, 
//Parallel, Anti-Parallel is derived from Image #3. An EDM 
would cause a polarity-dependent phase shift, A^, with 
the EDM given by d = h/S.(j)/ {AEt) (here r is the spin 
precession time with the E-field applied). An effect com¬ 
mon to both E-field polarities would produce an overall 
phase offset 0. Acj) was found to be uncorrelated with the 
other ht parameters. In both of the experimental runs, 
the measured EDM was found to be consistent with zero: 
— (4.0 ±5.2) X 10“^^ e-cm in the hrst measurement done 
with 3 mCi of ^^^Ra, and (0.6 ± 2.9) x 10“^^ e-cm in the 
second one, done with 6 mCi of ^^®Ra. The uncertainties 
listed above are statistical only. 

A variety of possible systematic effects could poten¬ 
tially lead to false EDM signals. Effects due to correla¬ 
tions between the high voltage and the probe frequency, 
the external B-held, the current supply for the bias B- 
held, and the standing-wave ODT power were consid¬ 
ered. Also modelled were those due to imperfect E-held 
reversal, induced B-helds due to the E-held pulsing, leak¬ 
age current, E x v effects. Stark interference in the ex¬ 
cited state hyperhne structure due to the presence of the 
AC laser helds and the DC E-held and geometric 
phase. The high voltage correlations were measured and, 
for each systematic effect, models were used to place an 
upper limit on the size of a potential false EDM. Using 
this analysis, all the systematic effects are calculated to 
be smaller than 1 x 10“^® e-cm. 

The effects due to imperfect E-held switching can also 
be constrained by direct measurements with trapped 
^^®Ra atoms. Effects quadratic in E-held do not produce 
an EDM-like signal, unless the E-held reversal is imper¬ 
fect. The voltage difference between the electrodes is 
recorded for each polarity; the E-held imbalance is found 
to be less than 0.7%. We use three E-held conhgura- 
tions (parallel, anti-parallel, and off) to directly measure 
both the linear and quadratic terms; the term that gives 
the phase a quadratic dependence on E-held allows us to 
place a limit on any potential E'^ systematics. The ability 
to place a limit on this effect is, therefore, limited by the 
statistical uncertainty of the spin precession ht — as this 
uncertainty improves, the limit placed on this effect will 
improve as well. For now, the quadratic effect was found 
to be below 2 x 10“^^ e-cm for the 3 mCi experimental 
run and below 5 x 10“^^ e-cm for the 6 mCi one. 

With the understanding that systematic effects are all 
negligible at the present level of sensitivity, we measure 
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FIG. 3. (Color online) Precession curves from the two experimental runs. Panels (a) and (b) are based on the first run, which 
used a 3 mCi ^^®Ra source, and panels (c) and (d) are based on the second run with 6 mCi. The panels show data with the 
E-field parallel to the B-field (red). E-field anti-parallel to the B-field (blue), and E-field off (black). Between the two runs 
the bias B-field was deliberately altered, resulting in the two different precession frequencies. An EDM would cause a phase 
difference between the E-field parallel and E-field anti-parallel fit curves. The global fits for the 3 mCi and 6 mCi runs yield 
X^/24 = 1.11 and X^/28 = 1.35, respectively. 


an EDM = -(0.5 ± 2.5stat ± 0.2syst) x IQ-^^ 

e-cm, and, hence, set an upper limit: |c?(^^®Ra)| < 5.0 x 
10“^^ e-cm (95% confidence). This represents the first 
EDM limit placed on an octupole-deformed species, and 
the first reported EDM measurement on atoms in an op¬ 
tical dipole trap. 

Improvements on the probing sensitivity by many or¬ 
ders of magnitude are possible. Detailed analysis indi¬ 
cates that the systematic effects of EDM measurements 
in an ODT can be reduced to below 10“^® e-cm [25] . 
Here, we discuss various approaches to improve statis¬ 
tical sensitivity. A factor of 10 in EDM sensitivity can 
be gained by increasing the standing-wave ODT lifetime, 
which can be achieved by improving the vacuum and by 
reducing heating. The gain in EDM sensitivity is linearly 
proportional to the gain in lifetime because the load- 
measurement cycle, currently 50 s, is dominated by the 
time needed for loading the MOT, and because an EDM 
measurement can be performed in parallel with MOT 
loading. 

The efficiency of detecting the spin state after the E- 
field has been turned off can be significantly improved 
by shelving one spin state to the long-lived ®Di level. 
Stimulated Raman Adiabatic Passage (STIRAP) [31] can 
be used to transfer atoms from the 7s^ ^Sq F = 1/2, 
mp = —1/2 state to 7s6d ®Di, while leaving those in the 
mp = -1-1/2 state undisturbed. This step turns the initial 
task of spin detection into atom-number detection with 
much higher signal-to-noise ratio, scattering 1000 pho¬ 
tons per atom instead of 3. At present, in shadow imag¬ 
ing, the detection noise is limited by the photon noise of 
the blue probe beam, independent of the atom number. 


Converting spin detection to atom-number measurement 
is expected to result in a factor-of-20 gain in EDM sen¬ 
sitivity. 

The trap can be made more efficient by using the 7s^ 
^So F = 1/2 —7s7p ^Pi F = 3/2 blue transition for trans¬ 
verse cooling and slowing. Under saturation, the laser 
force based on this transition is a factor of 100 stronger 
than the present scheme using the intercombination tran¬ 
sition. With this stronger laser force, trap loading effi¬ 
ciency can be improved by a factor of 30. This scheme, 
however, requires two additional repump lasers. 

Moreover, stronger sources of ^^®Ra are under devel¬ 
opment at various nuclear physics accelerator facilities, 
including the Facility for Rare Isotope Beams [32]. For 
example, it has been calculated that spallation of a tho¬ 
rium target induced by a 1 mA beam of deuterons at 1 
GeV will yield ^^®Ra at the rate of 10^® s“^ [33], which is 
5 orders of magnitude stronger than the currently avail¬ 
able supply. 

A ^^®Ra EDM experiment with an E-field of 100 
kV/cm, TV = 1 X 10® atoms, and r = 100 s of spin 
precession time can reach a statistical sensitivity at the 
level of 10~^ ® e-cm in T = 100 days, according to 
6d = h/(2Ey/TNT). Along the way towards this long¬ 
term goal, a recent analysis m suggests that a ^^®Ra 
EDM limit at the level of 10“^® e-cm would tighten the 
limits on certain types of T- and P-violating electron- 
nucleus interactions by at least an order of magnitude. 
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